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Abstract 



We show that the decaying magnetohydro dynamic turbulence leads to a 
more rapid growth of the correlation length of a primordial magnetic field than 
that caused by the expansion of the Universe. As an example, we consider 
the magnetic fields created during the electroweak phase transition. The 
expansion of the universe alone would yield a correlation length at the present 
epoch of 1 AU, whereas we find that the correlation length is likely of order 
100 AU, and cannot possibly be longer than 10 AU for non-helical fields. 
If the primordial field is strongly helical, the correlation length can be much 
larger, but we show that even in this case it cannot exceed 100 pc. All these 
estimates make it hard to believe that the observed galactic magnetic fields 
can result from the amplification of seed fields generated at the electroweak 
phase transition by the standard galactic dynamo. 
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I. INTRODUCTION 



Recently, considerable interest has been focused on the possibility that a primordial 
magnetic field may have been created at some early stage of the evolution of the Universe. 
While the existence of a weak widespread extragalactic chaotic magnetic field cannot be 
ruled out, almost the only place where such a primordial field may have left an observable 
imprint is in galaxies, many of which possess microgauss, kiloparsec scale fields that are 
thought to be the result of dynamo amplification of a weak seed field [|1|. The idea that the 
seed for the galactic dynamo may be the field created in the very early Universe has inspired 
a number of works. For example, such fields may have appeared at the electroweak phase 
transition ||. There are numerous other proposals involving physics at various scales (for a 
brief overview and further references, see ||). 

One feature shared by most particle-physics scenarios is the smallness of the correlation 
length of the magnetic field which results. Indeed, at the moment of creation, the correlation 
length is limited by the horizon radius. By pushing this moment to a very early stage in 
the history of the Universe, one makes the correlation length much smaller than it would 
be if the magnetic fields were created more recently, say, during proto-galactic contraction 
0]. For example, the fields generated at the electroweak phase transition, when the horizon 
radius is about 1 cm, would have a correlation length of at most about 10 15 cm in the present 
Universe. Contraction of proto-galaxies is likely to reduce this length scale by 2 orders of 
magnitude, which gives 10 13 cm, or 1 AU, as the characteristic length scale of the seed field 
which the galactic dynamo is supposed to amplify. Realistically, this length scale is likely 
much smaller, since the correlation length at creation is typically much less than the horizon 
sizeQ. 

Because any primordial magnetic fields generated at early times (via particle physics oc- 
curring at, say, the electroweak scale) have such small correlation lengths, these fields are not 
attractive candidates to play the role of seed fields for the galactic dynamo, unless the cor- 
relation length is somehow increased. The best developed theory of magnetic amplification 
in galaxies, the mean-field dynamo 0, operates under the assumption that the magnetic 
field is smooth on the scale of turbulent motion of the interstellar gas, which is of order 
100 pc. One may attempt to apply the mean-field dynamo theory for large-scale Fourier 
components of the chaotic magnetic field, neglecting the small-scale ones, however it is not 
obviously the correct procedure. That the magnetic field at very small scales exponentiates 
rapidly is a well known fact that poses a serious problem for the mean-field dynamo theory 
||. In the situation when the seed field itself resides at small scales, this problem is likely 
to become more severe. 

In this paper, we will not address the question of how a seed field is created, whether at 
the electroweak phase transition or during some other early epoch. Our goal is to investigate 



1 A magnetic field created during or before inflation || may have a large correlation length, however 
the proposed inflationary models that could possibly generate a large-scale magnetic field may seem 
rather contrived ||. We will not consider inflation-induced magnetic fields in this paper, but one 
should keep in mind this alternative. 
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the possibility that magnetohydrodynamic (MHD) effects can lead to a substantial increase 
in the length scale of the magnetic field at the present epoch. We show that decaying MHD 
turbulence typically leads to a faster growth of the magnetic correlation length than one 
would expect from the expansion of the Universe alone. We estimate that in the case of 
magnetic fields generated at the electroweak phase transition, the enhancement factor is 
10 2 , and the correlation length may reach 100 AU. If the primordial field has a large Chern- 
Simons number, the enhancement factor may be much larger, but the correlation length 
cannot possibly exceed 100 pc today if the magnetic field comes from the electroweak epoch. 
Thus, although the magnetohydrodynamic effects we consider are of some help, they cannot 
increase the correlation length enough to make the electroweak generation of a primordial 
seed field a viable option. We also consider generating the seed field at the QCD phase 
transition. This may be a possibility, but is only viable if the bubble separation at the phase 
transition is very large. 

This paper is organized as follows. In Sec. [TT] we write the basic equation governing 
the MHD of the Universe. Sec. |HI| is devoted to the decay of non-helical MHD turbulence, 
whereas Sec. [TV] describes the scaling laws of the decay of helical turbulence. Sec. [V] contains 
concluding remarks. The Appendix contains details about the EDQNM approximation used 
in our numerical simulation. 



II. THE MHD EQUATIONS 

To introduce notation, let us first write down the MHD equations for the early Universe. 
The MHD description applies since at interesting length scales the Universe is a good con- 
ductor. When the Universe expansion can be neglected (that is, when the evolution occurs 
on a time scale much smaller than the Hubble time), the MHD equations for a relativistic 
fluid have the form (for a derivation see, e.g., Rcf. ||), 

<9v 1 / Z? 2 \ 1 

+ ( v • V)v = — V [p + — + — — (B • V)B + z/V 2 v 

at P + P \ 2 J p + p 

— = V x (v x B) +r?V 2 B (1) 

where v is the velocity of bulk fluid motion, B is the magnetic field, p and p are the 
density and pressure of the fluid, respectively. We assume that the bulk flow velocity is non- 
relativistic, u < 1, although individual particles of the fluid move with the speed of light. 
The pressure can be eliminated from the equations by using the incompressibility condition, 
V ■ v = 00. The kinematic viscosity v and the resistivity rj, both having the dimension of 
length in units in which h = c = 1, are determined by microscopic physics. One can define 
the kinematic and magnetic Reynolds numbers, 



2 The fluid can be considered incompressible if the Mach number, i.e. the ratio between v and the 
sound speed, is much less than 1. Since for ultra-relativistic fluids the speed of sound is l/v3, this 
condition is essentially equivalent to v <C 1. 
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vl vl 
Re = — , Rgm = — 
v V 

where v and I are the typical velocity and typical length scale of the fluid motion. As a 
rule, in the early Universe the magnetic Reynolds number is much larger than the kinematic 
one (i.e. the magnetic Prandtl number, defined as v/r], is large.) Hydrodynamic turbulence 
occurs when Re ^> 1. For example, at T = 100 GeV, v ~ (aTlog^) -1 , and r\ ~ ^log^. 
Taking v ~ 10 _1 , I ~ 10 _2 a#, the two Reynolds numbers are Re ~ 10 11 and ReM ~ 10 17 , 
i.e. very large. 

Let us turn to the case of the expanding Universe. In the radiation-dominated epoch, 
neglecting the slow change of the number of degrees of freedom with temperature, the MHD 
equations are of almost the same form as Eqs. ([!]) , 

-V (p+^j+(B- V)B + z>Av 

V x (v x B) + r/AB (2) 

where t = Jdt a _1 (t) is the conformal time, B = B/^p + p is defined in such a way that the 
conservation of the magnetic flux during the expansion corresponds to the constancy of B. 
The spatial derivative in Eq. (0) V is taken with respect to the comoving coordinates. We 
have also defined v = a~ l v and fj = a _1 ?7, which are functions of time not only because of 
the factor a -1 but also due to the temperature (and hence time) dependence of the viscosity 
and conductivity of the Universe. The formal coincidence of the MHD equations in the 
radiation-dominated Universe and in the non-expanding Universe allows us to discuss the 
evolution of the magnetic field in both cases from a common point of view. 
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III. DECAY OF NON-HELICAL MHD TURBULENCE 

For definiteness, let us discuss the magnetic field generated during the electroweak phase 
transition, assuming the transition is first order [PjLOfl. It has been suggested in that after 
the phase transition, the cosmic fluid is in a state of turbulent motion, as parts of the fluid 
moving in front of different expanding bubbles collide p|. An important point that has not 
been sufficiently emphasized in the literature is that until either another first-order phase 
transition happens, or gravitational instability begins to act, no additional stirring takes 
place and the turbulence will decay freely by the turbulent energy dissipation (see below.) 
As a consequence, turbulence does not continue forever, but should terminate at some point. 
Therefore, to establish the growth law of the magnetic correlation length, one has to study 
the free decay of MHD turbulence, and its termination. 

Immediately after the phase transition, the characteristic velocity Vo of the turbulence 
is that of the bubble walls. The correlation length of fluid velocity is the bubble separation 
Iq. We will also assume the presence of an initial chaotic magnetic field that has the same 
correlation length as the velocity field and the same energy density as the turbulent bulk 
fluid motion. We rely here on the "equipartition" hypothesis, which states that any small 
seed magnetic field (which may have been created during the phase transition or come 
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from the hypercharge field created at some earlier epoch [fLOfl ) will be enhanced by turbulent 
convection until the magnetic energy is equal to the kinetic energy of the bulk fluid motion, 
Bq ~ Vq. Whether equipartition hypothesis is valid is a controversial issue. It is known 
that an arbitrarily weak magnetic field is exponentially enhanced by turbulence when ReM 
is larger than some critical value (of order 100, which is easily satisfied in cosmological 
conditions), but it is unclear what is the saturation level of the magnetic energy. However, 
it is very unlikely that we underestimate the magnetic field by relying on equipartition. As 
even this magnetic field will prove insufficient, we will not worry whether we are making an 
overestimate. (The values of vq and B mentioned above correspond to the integral scale, 
i.e. the largest scale of turbulence, which dominates the energy.) 

We will also presume for the present that the mean density of the Chern-Simons number 
is vanishing. This is the case if the initial condition is parity symmetric in the statistical 
sense. We discuss the case of nonzero Chern-Simons number, and why it is different, in a 
later section of this paper. 



A. Non-expanding Universe: scaling laws 

Let us first consider the simpler case, in which the expansion of the Universe can be 
neglected. This corresponds to the first step of the evolution, when the time passed from 
the creation of the magnetic field is smaller than the Hubble time. We will first make 
a simple argument of how the correlation length grows with time, and then discuss the 
assumptions behind this argument. The mechanism for the growth is selective decay: modes 
with larger wavenumbers decay faster than those whose wavenumbers are small. Thus, at 
increasing time only modes with smaller and smaller wavenumbers survive. Let us assume 
that at some time t only modes with length scale larger than / survive. The modes that 
are about to decay, i.e. those at the scale I, dominate the kinetic and magnetic energies. 
Denote the typical value of v and B of these modes as v\ and B\. From the equipartition 
hypothesis, one has Bi ~ vi. Let us find the dependence of v\ and Bi on I. If we assume that 
modes with wavelengths larger than I remain the same during the decay of shorter modes, 
the configuration of v and B at moment t is essentially the initial configuration, smeared 
out over the length scale /. At the time moment t, the value of v and B at a point in space 
is roughly the average of the initial v and B over the sphere of radius I around this point. 
Assuming that v and B, at t = 0, is not correlated over lengths larger than Iq (which is 
certainly true if l is the Hubble size of the Universe at this time), the average is over (///o) 3 
patches with random v and B. Therefore, the order of magnitude of v and B at the moment 
t is 

The characteristic time for the decay of modes with length scale I is the eddy turnover time, 

n ~ - (4) 

Vl 

Notice that the decay time (f|) does not depend on any viscosity and resistivity. This is the 
consequence of the turbulent cascade: during the time n the energy is transfered from modes 
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with length scale I to 1/2, then it is cascaded to 1/4, 1/8 etc. with faster and faster rate until 
the viscosities enter the scene at some small scale Zdiss (which is called the dissipation scale). 
Due to the cascade, the energy (both kinetic and magnetic) is dissipated at a much faster 
rate than it would be without turbulence. The length scale /di SS does not enter any of our 
formulae, as we are interested in the much longer length scales at which the cascade begins; 
all that matters is that dissipation occurs at some very small length scale, terminating the 
cascade. 

Since I is the smallest wavelength that survives at the moment t, its lifetime must be 
comparable to t. Solving 77 ~ t using Eq. (^), one finds that the correlation length scales 
with time as, 

/ t \ 2/5 

l ° ( 5 ) 



to J 

where t = Iq/vq is the eddy turnover time at t = 0. From Eqs. (QJ|) and the fact that the 
energy is proportional to v 2 and B 2 , we find that the latter decays as £~ 6 / 5 . It is worth noting 
that all our arguments leading to these scaling laws are also valid for decaying hydrodynamic 
(not MHD) turbulence. Indeed, the t _6//5 decay of the turbulent energy is known in fluid 



dynamics as Saffman's law [pj , p^| . 



The implicit assumption behind our derivation of the t 2//5 law (^j) is that the spectrum 
at low k remains unmodified during the decay of high-fc modes. Such an assumption is 
essential for deriving Eq. (^) (called "volume averaging" elsewhere.) The validity of this 
assumption may seem uncertain, since inverse cascade, i.e. the process of transferring energy 
from high to low modes happens under some conditions in MHD. In |l3|,p|, it was suggested 
that inverse cascade is a generic feature of MHD turbulence and could greatly enhance the 
correlation length of the magnetic field (although, coincident ally, the scaling law for the 
growth of the correlation function obtained in is the same as our Eq. (|5|).) What 

we want to argue here, to support the derived scaling laws, is that the inverse cascade 
of magnetic energy required a special condition to be met, namely, the average density of 
the Chern-Simons number must be non-zero. That the conditions for the inverse cascade 
are rather restrictive is a well known fact in plasma physics: previous studies fl4Hl6f have 
identified these conditions to be either a non-zero Chern-Simons number density of the 
magnetic field, or a non- helical turbulent flow (see also the discussion below.) Since we do 
not expect the bulk fluid flow to be helical, the only chance to have an inverse cascade is if 
the magnetic field has non-zero Chern-Simons number, which is not the case considered in 
this Section. We have also checked that there is indeed no inverse cascade in the non-helical 
MHD turbulence decay by simulating it using the eddy-damped quasi-normal Markovian 
(EDQNM) approximation. Our simulation will be described in the next subsection. 



B. Numerical simulation of decaying non- helical turbulence in a non-expanding 

Universe 

To verify the absence of the inverse cascade, as well as to check our scaling laws, we use 
the EDQNM approximation for numerical simualations. A full discussion of the merits and 
shortcomings of the EDQNM approximation can be found in Refs. |T2|,|T3, here we will only 
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note that this approximation is a member in a family of closure schemes, which are similar 
to the Hartree-Fock approximation in field theory. The EDQNM approximation has been 
found very successful in reproducing qualitative features of hydrodynamics and MHD as 
the Kolmogorov spectrum of driven hydrodynamic turbulence, the inverse cascade in MHD 
(under the above-mentioned conditions), etc. One should be aware that it is unclear how the 
EDQNM approximation can be systematically improved. However, it should be adequate 
to answer our qualitative questions. 

In the closure schemes, the turbulence is characterized by two functions E k and M k , 
which are the spectra of kinetic and magnetic energy (in the case of helical flow and helical 
magnetic fields, one should introduce two more functions of k.) The full set of EDQNM 
equations are quite complicated In the Appendix, we reproduce the equations for the 
case of non-helical turbulence. As the initial condition at t = 0, we choose the initial kinetic 
and magnetic spectra to have the same functional form, 



The k 2 behavior of E k and at low k's comes from the assumption that the initial state 
consists of uncorrelated eddies^. In the initial condition (|6|), the magnetic energy constitutes 
10% of the initial total energy at t — 0, which is equal to 1. In choosing the initial magnetic 
energy to be small compare to the total energy, we obtain as a bonus the possibility of 
checking whether equipartition is reached during the decaying turbulence. We choose v = 
7] = 10~ 3 , which correspond to the Reynolds numbers of order 10 3 at t = 0. This is much 
smaller than the cosmological Reynolds numbers, but large enough for the physics at the 
integral scale to be insensitive of the exact values of Re and ReM- We use logarithmic 
discretization in the k space, with 4 points per octave (we have checked that doubling the 
number of points per octave does not alter the result in any substantial way), and the 
simplest first-order time-update algorithm. 

The main results are summarized in the Fig. ([T|-§|). In Fig. [I], the spectra of kinematic 
and magnetic degrees of freedom are plotted at t — and t = 100 (measured in units of the 
eddy turnover time at t — 0). We see that the low-/c part of the spectrum is not modified by 
turbulent decay, which means that there is no inverse energy cascade. It is interesting to note 
that, at if: = 100, in the high-fc part of the spectrum the magnetic energy exceeds the kinetic 
energy, and there is a range of k where the Iroshnikov-Kraichnan spectrum Ef* ~ M k ~ k~ 3 / 2 
is observed. Both features have been seen before in the simulation of forced (not decaying) 



3 To see that, one could use the following argument. The typical wavelength in the initial state is 
ko ^ Iq 1 - If one throws away all modes with momentum larger than some k <C ko, this is equivalent 
to smearing out the velocity and the magnetic field over the distance of k^ . As discussed before, 
if v and B do not correlate over lengths larger than k$ , the average takes both down by a factor of 
(/c//co) 3//2 . The energy, which is proportional to v 2 and B 2 , is down by a factor of (k/ko) 3 . Compare 
this with the case of an energy spectrum with k 2 behavior at low k, in which throwing away modes 
with energy larger than k makes the energy of remaining modes proportional to f^dk k 2 ~ k 3 . 



E k = 0.9- 




M k = 0.1- 




(6) 
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turbulence |1| . In Fig. 0, the typical wavelengths of kinetic and magnetic modes, which are 
defined as 



fdpE p JdpM p 
Jdpp^Ep 1 jdpp~ l M p 

are plotted versus time from t — to 1000. Both k(t) and ku(t) decreases at large t and 
their time dependence agrees very well with our prediction, k ~ t -2 / 5 . The ratio k^(t)/k(t) 
approaches a finite value (about 2) at large t. We conclude that our simple arguments 
reproduce correctly the behavior of the correlation length in the non-helical turbulence decay. 
In Fig. H we plot the dependence of the total kinetic and magnetic energy as functions of 
time. The long-time behavior is well described by Saffman's law E ~ t _6//5 . Notice that 
the magnetic energy grows for a few eddy turnover periods before decaying. Apparently, 
equipartition is at work here: at first some energy is pumped from the kinetic motion into the 
magnetic field. Equipartition is never fully realized, but the ratio of the magnetic energy to 
the kinetic energy approaches a value closed to 1 (namely, 0.8) in our simulation. Therefore, 
we have seen that our qualitative arguments (the most important being the lack of inverse 
cascade) agree well with the numerical simulation. 



C. Expanding Universe 

So far we have considered the case of a non-expanding Universe. In the radiation- 
dominated Universe, the correlation length, in comoving coordinates, grows as r 2//5 , where 
t is the conformal time. Since r ~ T _1 (T being the temperature), we conclude that in 
this case the growth of the correlation length follows the T~ 2 / 5 law (not counting the trivial 
redshift). 

In the real situation we have a hybrid of the two cases. If at t = the eddy turnover 
time to ~ lo/vo is smaller than the age of the Universe, there will be a substantial period 
of time when the expansion of the Universe can be neglected, that is, when < t ^ tu- 
The expansion factor for the length is then (tHVo/lo) 2 ^ 5 . When t ^ tn, one has to take 
into account the expansion, and the increase in the correlation length (not counting the 
expansion) is (Tew/To) 2//5 , where To is the temperature of the Universe when the fluid flow 
becomes non-turbulent (we determine To below.) Therefore, the final correlation length at 
the present epoch is, 

, , ft H v \ 2/5 /T EW \ 2/5 T EW 

fnow = to • I —j I • I -7f7~ ) ■ 7f , {') 

\ to / \ J-0 / J- now 

where the last factor comes from the expansion of the Universe. Putting in numbers, using 
for T the temperature of e + e~ annihilation (which is of order m e ), using Vq ~ 0.1 and 
lo ~ 10~ 2 a#, one finds Z now ~ 100 AU. This scale is larger than the scale one would have 
if the only mechanism of dissipation of magnetic energy is resistive diffusion (about 1 AU 
pTp, but it is still too small to be a direct seed for the mean-field dynamo. The situation 
is better if one takes larger values of Vq and lo. However, even in the extreme case when 
the magnetic field is initially correlated on the horizon size and its initial energy density is 
comparable to that of the whole Universe, the correlation length is still less than 10 4 AU. 
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Let us check that To is the temperature at which e + e~ annihilation occurs, which is to 
say that below this temperature, turbulence ceases, and therefore so does the enhancement 
of the magnetic correlation length. To this end, we need to compute the Reynolds number 
Re before and after the e + e~ annihilation. At T e + e - ~ 1 MeV, using Eq. ([?]) and replacing 
T now by T e + e -, one finds the correlation length I ~ 10 19 GeV -1 . The velocity at T is 

V ~ Vq • — — • — — ~ 10 



and the viscosity is v ~ (a 2 T) _1 ~ 10 7 GeV _1 . The kinematic Reynolds number is then 10 5 . 

After e + e~ annihilation, however, the cosmic fluid becomes much more viscous. This is 
the result of the dramatic drop of the density of charged particles, which leads to a very 
long mean free path for photons, which now dominate the viscosity. The viscosity jumps by 
the photon-to-baryon ratio, 10 10 , and the Reynolds number becomes small. Therefore, e + e~ 
annihilation is the moment when turbulence terminates^. After this moment, the magnetic 
field is frozen into the essentially non-moving fluid and the correlation length increases only 
by the expansion of the Universe, as we assumed in deriving Eq. (0). 



IV. DECAY OF HELICAL TURBULENCE 

Let us now turn to the case when the original magnetic field is helical. Here helicity 
refers to the Abelian Chern-Simons number, Acs — Jdx. A • B, which is called the magnetic 
helicity in plasma physics. The situation when the initial field has a non-vanishing Chern- 
Simons number is not entirely hypothetical, since some mechanisms indeed generate such 
helical fields [[RJ . However, our interest will be mostly in finding out the largest length scale 
one could possibly get from a magnetic field from the electroweak epoch. 

It is a well known fact in plasma physics that the Chern-Simons number is an approx- 



imately conserved quantity of MHD turbulence |L4| . If the fluid is an ideal conductor, the 
conservation of Acs is exact. This can be seen from the geometrical interpretation of the 
Chern-Simons number: Acs is proportional to the number of links of the magnetic field 
lines. When rj = 0, the latter are frozen in the moving fluid, and the number of links cannot 
change with time. At nonzero rj, the conservation of Acs is less obvious, since another 
conserved quantity, the energy E = |/<ix((p + p)v 2 + B 2 ), dissipates at a finite rate at 
arbitrarily small v and rj. The smallness of the dissipation of Acs can be explained in the 
following (non-rigorous) way, which relies only on the difference between the dimensions 
of E and Acs- Let us write down the equations for the dissipation of the energy and the 
Chern-Simons number, which can be derived from Eq. (|]) fl4| , 



4 It is sometimes, based on the results of [18|, argued that the large diffusion length of photons leads 
to damping of the magnetic field, similar to the Silk damping of density fluctuations. However, 
Ref. [18] deals with the damping of MHD waves but not the background magnetic field on top of 
which these waves propagate (see also the analysis in Ref. fl9|l.) 
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E = — Jdx (v(p + p)(V x v) 2 + x B) 2 ) (8) 
N C s = — 2?7 fdxB'(VxB) (9) 



Let us compare the order of magnitude of the right hand sides of Eqs. (HH), assuming 
for simplicity v ~ r\ and equipartition, pv 2 ~ -B 2 . If both i£ and Acs are dominated by the 
largest length scale of turbulence (the integral scale) /, then Nqs/E ~ /. On the other hand, 
the dissipation of both energy and the Chern-Simons number occurs at the dissipation scale 
/diss- Since /diss is much smaller than / (by a power of the Reynolds numbers), we conclude 
that Nqs/E <C Nqs/E, which means that the dissipation rate of the Chern-Simons number 
goes to zero in the limit of small viscosity and conductivity, in contrast to the energy which 
dissipates with a finite rate in this limit. When the Reynolds numbers are large, one can 
neglect the dissipation of the Chern-Simons number (but not of the energy.) 

The conservation of Acs has an important consequence for the evolution of the magnetic 
field. When Acs is non-vanishing, the short-scale modes are not simply washed out during 
the decay: their magnetic helicity must be transfered to the long-scale ones. Along with the 
magnetic helicity, some magnetic energy is also saved from turbulent decay. This process 
is known as the inverse cascade [14,15][]. The coupling with kinetic motion ensures that a 
part of the energy cascaded to large scales also goes to the kinetic energy of fluid motion, 
maintaining equipartition. We will assume maximal helicity, i.e. A • B has the same sign in 
all space, so the magnetic helicity density is proportional to AB ~ IB 2 . Since this quantity 
is conserved during the decay of the turbulence, the field strength should scale with the 
correlation length as 

B { - 

This law replaces Eq. (0). This corresponds to "line averaging" that gives a much larger 
amplitude of the magnetic field than the "volume averaging pi . However, our averaging is 




not postulated, but is a dynamical effect which occurs only in the helical case. An argument 
similar to the one used in the non-helical case now gives the following decay law for helical 
MHD turbulence: 

/ t \ 2/3 
' ~ *o ( - ) , 



t J ' 

which replaces the formula (||). Therefore, the correlation length grows faster if the initial 
field configuration is helical. Obviously, this comes from a slower decay of the magnetic field, 
which implies a larger equipartition velocity and hence a shorter relaxation time compared 
to the non-helical case. 



5 Inverse cascade occurs also in two-dimensional hydrodynamics and MHD. In both cases, the con- 
served quantities governing the inverse cascade, the enstrophy and the mean square magnetic po- 
tential 1 14], are positively defined, so the inverse cascade does not requires special initial conditions. 
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To see what is the largest magnetic length one can get, let us assume that at T = 100 
GeV one has l ~ a H and v ~ 1. Assuming the turbulence has all decayed at T , the 
correlation length at T = T is 



I = lr 



Tew \ 5/3 
T 



The velocity now scales as v ~ (Tew/To)- 1 / 3 . In the helical case, with our choice of parame- 
ters, the turbulence in fact survives beyond e + e~ annihilation. When T -C m e , the viscosity 
of the Universe is mostly due to photons and is given by 
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n e a T r] B aT 3 



where n e is the electron density, t\b is the baryon-to-photon ratio, and <jt ~ am~ 2 is the 
Thompson cross section. One sees that the kinematic Reynolds number drops below 1 at 
T ~ 100 eV. Therefore, the correlation length is now 

/T EW \ 5 / 3 / T \ 
I ~ l • -777- ■ - — ~ 100 pc 

Since the collapse of proto-galaxies will reduce this correlation length by a factor of 100, 
the galactic seed field will be correlated on a scale of order 1 pc, still much smaller than 
the scale of the turbulent motion of the interstellar gas which we recall is of order 100 
pc. Therefore, the presence of helicity in the initial magnetic field strongly enhances the 
magnetic correlation length, but still cannot make it large enough to be the seed for the 
mean-field galactic dynamo. 



V. CONCLUSION 

In this paper we have considered the behavior of the correlation length of the primordial 
magnetic field, if such a field is generated in the early Universe. We found the decay law 
of the MHD turbulence, which is responsible for a substantial increase of the correlation 
length. We also observed a qualitative difference between the cases of non-helical and helical 
initial magnetic field. Taking the case in which the magnetic field is generated during the 
electro weak phase transition as an example, we find that if the magnetic field is not helical 
the factor gained from MHD is about 10 2 and today the field may be correlated at scales 
as large as 100 AU. If by some chance the field is helical, the enhancement factor is much 
larger, but the final correlation length cannot exceed 100 pc, corresponding to 1 pc after 
proto-galactic collapse. 

Let us note some uncertainties remaining in our estimation. (In general, addressing these 
issues would lead to less optimistic values of the enhancement factor.) First, it is not at all 
clear whether in MHD turbulence the magnetic energy reaches equipartition with the kinetic 
energy of bulk fluid motion. While semi-analytical calculations (including our simulation) 
favor equipartition | [15|] , some numerical results indicate that the mean magnetic energy 



density remains small and concentrated at scales shorter than the largest turbulence scale 
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T6| . If the latter remains true at very high Reynolds number, the magnetic field will be 



smaller after the electroweak phase transition and the correlation length at present will be 
smaller than we have estimated. 

The second factor is neutrino diffusion. After the electroweak phase transition, neutrinos 
are the particles with the longest mean free path. There is a certain time interval when 
neutrinos are still not decoupled from the physics at the turbulence scale, but their diffusion 
length is so large that the neutrino contribution to viscosity makes the Reynolds number to 
drop below 1. During this time, there is no turbulence and the magnetic field is frozen. After 
the neutrinos decouple from the fluid motion at the turbulence scale, magnetic stress leads 
to restoration of turbulence, and the magnetic length may continue to grow. The overall 
effect is some reduction in the estimate of the final magnetic correlation length. 

It would be nice to end this paper on a more positive note, and to this end we turn 
now to investigating the QCD phase transition, which occurs later than the electroweak 
transition and so may yield a longer correlation length. If the QCD phase transition is first 
order, it will introduce fresh turbulence stronger than the decaying electroweak turbulence. 
The amplitude and correlation length of the magnetic field will be determined by the QCD 
turbulence. If fields were regenerated at the QCD transition on length scales of order of the 
horizon, the length scale it would have today is very large, since the horizon size at the QCD 
phase transition is much larger than at the electroweak epoch. Even with no enhancement 
from MHD effects, the expansion of the universe yields a magnetic correlation length of 
order 1 pc. It can be estimated that the turbulence will survive till at least matter-radiation 
decoupling, at which the correlation length is is of order 1 kpc in the non-helical case and 
100 kpc if the field is helical. However, a magnetic field correlated on the horizon size is 
not expected to be produced at the QCD phase transition, and these estimates are too 
optimistic. Instead, the initial length scale must be of the order of the bubble spacing at the 
end of the phase transition (which is the natural scale of turbulence). Bubble spacings larger 
than 100 cm are unlikely PI[ and may have problem with standard big bang nucleosynthesis 



22]. A bubble spacing smaller than 100 cm again implies a very small correlation length 



at the present epoch. If this constraint is respected, our analysis of the QCD case must 
end pessimistically, as in the electroweak case. However, one cannot completely rule out the 
possibility of a very large bubble spacing, which leads to a magnetic field correlated on a 
long enough length scale to be of interest. A nonstandard evolution in which mixing after 
the phase transition occurs rapidly by hydrodynamic flows instead of slowly by diffusion or 
nonstandard nucleosynthesis may be required in this case. 

Finally, although in this paper we presented the smallness of the magnetic correlation 
length as undesirable and tried to overcome it by invoking MHD turbulence, there may 
exist a non-standard dynamo mechanism where a small-scale seed gives rise to a large-scale 
magnetic field (see, e.g., f23|). This possibility is, however, outside the scope of the present 
paper. 
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APPENDIX A: THE EDQNM APPROXIMATION 

As mentioned in the main text, in closure schemes, the non- helical turbulence is char- 
acterized by the kinetic and magnetic energy spectra, E k and M k , defined so that the total 
kinetic and magnetic energies are J£°dkE k and J£°dkM k , respectively (we assume isotropy, 
so the spectra depend only on the absolute value of k.) The evolution of the spectra, in the 
EDQNM approximation, are described by the following equations, 



(d t + 2uk 2 )E k = Jdpdq6 kpq 



-(xy + z 3 )(k 2 E p E q - p 2 E q E k ) - V -z{\ - y 2 )M q E k + 



Q 



(d t + 2 V k 2 )M k 



+—z(l - y 2 )M„M ( 
q 

[k 3 

dp dq 9 kpq 



kp 



1 + xyz)M p E q - -qi - y 2 )E q M k - 
pq q 



k 2 



--z(l-y 2 )M q M k 



(Al) 



which are the non- helical version of the full set of EDQNM equations that can be found in 
Ref. ||15|| . In Eq. (|A"1|), the integration over p and q is perform in the region where k, p and 
q form a triangle; x, y and z are the cosines of the angles of this triangle, 



x 



p 2 + q 2 — k 2 



2pq 



q 2 + k 2 — p 2 
2qk 



k 2 + p 2 — q 2 
2kp~ 



and 9 kpq is defined as 



'kpq 



f^kpq 



where 



l^kpq j-^k l^p f^q 



[v + rj)k 2 + C s (£dqq 2 (E q + M q ]j' 2 + C a k(J*dqM K 



f^k 



1/2 



C s and C a are phenomenological constants, the exact values of which are not qualitatively 
important. Following we choose C s = 0.26 and C a = l/y/3. 
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FIG. 1. The kinetic and magnetic spectra at t = and t = 100 
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FIG. 2. Time evolution of the typical wavelength 
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